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Abstract 

We theoretically investigate the effects of laser polarization on photoelectron angular distribution 
through laser-induced continuum structure. We focus on a polarization geometry where the probe 
and dressing lasers are both linearly polarized, and change the relative polarization angle between 
them. We find that the total ionization yield and the branching ratio into different ionization 
channels change as a function of the relative polarization angle, and accordingly the photoelectron 
angular distribution is altered. We present specific results for the ^Pi/2-®Pi/2 an d 4j?3/2~6p3/2 
systems of the K atom, and show that the change of the polarization angle leads to the significant 
modification of photoelectron angular distribution. 
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I. INTRODUCTION 



The similarity between the autoionizing structure (AIS) and the laser induced continuum 
structure (LICS) is very well known. In an AIS process a discrete state lying above the ion- 
ization threshold is coupled to the continuum through configuration interaction, a resonance 
structure being created. In LICS two bound states are coupled to the common continuum 
through two laser fields (the probe and dressing lasers). The bound state coupled to the 
continuum through a strong laser (dressing laser) can also create a resonance structure hav- 
ing AlS-like properties, but compared to the AIS resonance its position and width are now 
controllable by the frequencies and intensities of the lasers. The first experimental observa- 
tion of LICS was successfully reported in Refs. . More comprehensive information on 
LICS can be found in a review paper by Knight et al 

Through LICS not only the ionization yield but also a number of some other processes 
can be altered: Several works based on LICS investigated non-linear optical effects such 
as the enhancement of third-high harmonic generation Q. In Ref. Q, effects of LICS on 
spin-polarization were studied for heavy alkali atoms. Recently, LICS for multiple continua 
was experimentally and theoretically investigated fl. In Ref. fl, the control of ionization 
products in LICS was suggested for the case of decay into multiple continua. 

It is well known that the photoelectron angular distribution (PAD) provides more infor- 
mation about the ionization process than the angle-integrated ionization signal [8]. PAD's 
of Na by the two linearly polarized lasers with a variable polarization were reported in Ref. 
{9], and the phase difference between the continua with same parity was extracted. By 
measuring PAD's of an alkali atom in a bichromatic laser field, a theoretical method was 
proposed in Ref. [lOj in order to extract the phase difference of the continua with opposite 



parities. 

n 

Most recently we have theoretically investigated how LICS affects PAD and specific 
results have been presented for the K atom. In Ref. however, we have assumed that the 
probe and dressing lasers are linearly polarized along the same direction. A natural question 
would be how PAD is modified, through LICS, by changing the relative polarization angle 
between the probe and dressing lasers. 

I — I 

The aim of the present paper is to generalize our previous work We now vary the 
relative polarization angle and analyze the modifications of the ionization yield and the pho- 
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toelectron angular distribution through LICS, and see how the enhancement or suppression 
of a particular ionization channel takes place. Since the photoelectrons ejected into different 
involved continua have different angular distributions, and those angular distributions de- 
pend on the relative polarization angle, we expect important modifications in terms of the 
ionization yields, branching ratios, and photoelectron angular distributions. 

The paper is organized as follows. In Sec. II we present the theoretical model: The 
time-dependent amplitude equations which describe the dynamics of the LICS process are 
derived, and the ionization yield and photoelectron angular distribution are calculated. The 
consistency of our results has been checked using an alternative approach based on the 
density matrix equations. The theoretical results obtained using these two formalisms are 
of course identical. Further details about the density matrix formalism are provided at the 
beginning of Sec. Ill, and in Appendix |A] Sec. Ill is mainly devoted to the discussions on 
the numerical results for the total and partial ionization rates, branching ratios, and PAD's 
at different polarization angles. The atomic parameters needed for the LICS calculation are 
given in Appendices iBllDl 

II. THEORY 

The system we consider in this paper consists of an initially occupied 4p state, initially 
unoccupied Qp state, and the continuum states of the K atom together with the linearly 
polarized probe and dressing lasers that couple 4p and 6p to the continuum states, respec- 
tively. This implies that, prior to the interaction of the system with the probe and dressing 
lasers, K atoms in the ground 4s state have been excited to the 4p state by a linearly po- 
larized auxiliary laser. By choosing an appropriate frequency of the auxiliary laser, we can 
selectively excite either Api/2 (rrij = ±1/2) or 4^3/2 (rrij = ±1/2), which will serve as an 
initial state in this work. For simplicity, we assume that the polarization axis of the dressing 
laser is parallel to that of the auxiliary laser, while the polarization axis of the probe laser 
can be arbitrary. 

Here we are interested in a particular geometry where the polarization vector of the 
dressing and auxiliary lasers are along the z— axis and that of the probe laser is assumed 
to lie in the xz-plane, as shown in Fig. [U For such a case, the polarization vector of the 
probe laser is defined as = eisin^p + e3Cos# p , where ei, and e3 are the unit vectors 
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along the x and z axes, respectively, and 9 P represents the polarization angle of the probe 
laser with respect to that of the dressing laser. Defining the frequencies of the probe and 
dressing lasers as co p and Ud, respectively, the total electric field vector can be written as, 

E(i)= £ £ 0a {t)e^cos{u a t). (1) 

a=p,d 

A Gaussian temporal profile was employed for the amplitude of the laser fields: £o a (t) = 
So a exp —4 In 2 (t/r a ) 2 , where r a represents the temporal width for the full width at half 
maximum (FWHM) of the probe or dressing pulse with a = p or d, indicating the probe 
and the dressing pulses, is the polarization vector of the laser pulse a. 

Based on the above assumptions, the level scheme we consider in this paper is now 
described in Figs. El(a)-I2(c) for the K 4pi/ 2 -6pi/2 system, at 9 P — 0°, 9 P — 90°, and in 
between, i.e., 0° < 9 P < 90°. If both polarization axes of the probe and dressing lasers are 
parallel, i.e. 9 P = 0° as shown in Fig. [2(a), due to the selection rule my = rrij (where the 
prime index is used for quantum numbers of the continuum) the entire 4pi/2-6pi/2 system 
with rrij = ±1/2 can be decomposed into the two independent subsystems which consist of 
4pi/2 ( m j = ±1/2), 6pi/2 {rrij = ±1/2), and the continua es {my = ±1/2) and ed {my = 
±1/2). The ionization yields for these subsystems are obviously symmetric to each other, 
and for simplicity we can consider only one of them, as already explained in our previous 
paper llj. Similarly, at 9 P = 90°, because of the selection rule my = mj ± 1, the entire 
4pi/2-6pi/2 system with mj = ±1/2 can be decomposed into the two independent subsystems 
consisting of ^pi/ 2 {mj = ±1/2), 6pi/ 2 {mj = =(=1/2), and the continua es (m^/ = ±1/2) and 
ed {my = ±1/2) with additional incoherent channels ed {my = ±3/2), as shown in Fig. 
H(b). Again, both subsystems are completely symmetric, and it is sufficient to study only 
one of the two subsystems. However, for the intermediate values of the polarization angle, 
0° < 9 p < 90° as shown in Fig. [2(c), because of the selection rules my = mj ± 1 (from the 
perpendicular component of the probe polarization vector with respect to the quantization 
axis, 9 P = 90°) and my = mj (from the parallel component, 9 P = 0°), the entire system 
4pi/2-6pi/2 cannot be decomposed into the two independent subsystems anymore, and the 
entire system, 4p!/ 2 , Qpi/2, and the continua with all possible magnetic sub levels, have to be 
taken into account at the same time. The continuum states |q,), {b = 5,8), not presented 
in Fig. [2J have the same quantum numbers as the continuum states \c a ), {a = 1,4, for 
my = ±3/2, ±1/2, —1/2, —3/2), but they correspond to a different value of energy because 



of the incoherent one-photon ionization from 6pi/ 2 by the probe laser. A similar level scheme, 
taking into account appropriate dipole selection rules, has been considered for the K 4p 3 / 2 - 
6j>3/2 system. 

In order to observe LICS, it is necessary that the initially occupied 4pi/ 2 (or 4p 3 / 2 ) state 
(denoted as |1) for the magnetic sublevel having rrij = +1/2 and |2) for the magnetic sublevel 
having rrij = —1/2), and the initially unoccupied Qpi/2 (or 6^3/2) state (denoted as |3) for 
rrij = +1/2 and |4) for rrij = —1/2) are coupled by the probe and dressing lasers whose 
frequencies nearly satisfy the two-photon resonance, i.e., E ip + u p ~ E ep + Ud- As long as we 
use a ns laser with appropriate intensities and detunings, it is perfectly valid to treat each 
4pi/2 — Gpi/2 and 4p 3 / 2 — 6p 3 / 2 system separately, as it was explained in our previous paper 



It should be mentioned that we expect a different behavior of the two systems: For 
the K 4p 1 / 2 — 6pi/2 system the initial state, 4pi/ 2 , is an isotropic mixture of all possible 
magnetic sub levels (recall that the magnetic sublevels rrij = ±1/2 are equally populated by 
the auxiliary laser), implying that the initial state is spherically symmetric. It is obvious that 
the PAD from the spherically symmetric initial state orientates along the polarization axis 
of the probe laser if the dressing laser is off. PAD changes neither its shape nor magnitude 
12|. That is not the case for the K 4p 3 /2 — 6p 3 /2 system, since not all possible magnetic 
sublevels are excited with the same probability, and accordingly the initial state, 4p 3 / 2 , is 
non-spherical (polarized). Therefore, we expect quite different modification of PAD for the 
4pi/2 — Gpi/2 and 4p 3 / 2 — 6p 3 /2 systems through LICS, by varying the relative polarization 

n 

angle, as we have already seen for 9 P = 0° [ll|. Further details will be provided in Sec. III. 
A. Time-dependent amplitude equations 

In order to study the temporal evolution of the atomic system in laser field, we have 
used the standard procedure as described in our previous paper [TJJ. Briefly, we solve the 
following set of time-dependent amplitude equations: 

«i = -^f 1U1 - itti 3 (1 - — ] u 3 - iQ 14 , 1 1 - — ] « 4 , (2) 

2 V ^13/ V 914/ 

ii2 = --f 2 u 2 - *n 23 (l — ) u 3 - iQ 2 A (l — ) u 4 , (3) 

2 V 923/ V 924 / 
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u 3 = 



W4 = 




(5) 



(4) 



where u/s (j = 1,4) represent the probability amplitudes of states Note that all the 
probability amplitudes for the continuum states have already been adiabatically eliminated 
in the Eqs. (]2j)- (JHJ) . 5 is a two-photon detuning defined by 5 = S sta uc + S stark, where the static 
detuning is defined by 8 static = (Ei + hu> p ) — (E 3 + hud), and S sta rk is a total dynamic ac Stark 
shift defined by 5 sta rk — (S{ P ^ + S[ ) — (S 3 P ^ + S 3 d ^). In all the numerical results presented in 
this work the zero point of the detuning has been chosen such that 5 — > 5 — S^l^ k , since the 
ac Stark shifts simply translate the LICS resonance on the detuning scale. The superscript 
of S™%rk means that the ac Stark shift is calculated at the peak value of laser intensity D^) 's 
are the bound- free matrix elements by laser a (a = p or d) from the bound state \j) to the 
continuum |c), which are connected to the partial ionization widths through the relation 
r^f = 2ir\Dj^\ 2 . Tj represents the total ionization width of state \j), i.e., Tj = 7j + Tj 
(for j = 1,2) and = jj + 1^ + 1^ (for j = 3,4), where 7j is the phenomenologically 
introduced spontaneous decay width of state In the above equations the two-photon 
Rabi frequency, flij, can be written as a sum of the partial two-photon Rabi frequencies into 
the coherent es and ed continua of energy e, i.e., 



where and represent the total and partial asymmetry parameters, respectively. It is 
very well known that, for the light alkali-metals such as Li, N, and K, the dependence of 
radial matrix elements and phase shifts on the total angular momentum quantum number 
j is very small and it can be neglected [3]. Now, the following relations are satisfied by 



the angle-integrated atomic parameters: 1^ = , I3 = I4 , S[ aS) = Sz \ S 3 a ^ = S/?\ 



^13 = ^24, ^14 = — ^23, <?i3 = <?24, and gi 4 = —^23- Details about the calculation of 
the atomic parameters such as Rabi frequencies, ionization widths, ac Stark shifts and 
asymmetry parameters are given in the Appendices [BllDl 

Since the behavior of the population dynamics in the continuum is of our interest, we 
also need the following set of amplitude equations for the continua: 





u, 



-ticaUca - 




2 U 2 - 




3 U 3 - lLJ Ca 



4 u i , 



(7) 
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u 



i8 Cb u Cb - iD ( £lu 3 - iD%{u 4 



(8) 



Here u Ca represents the probability amplitude of the coherent continuum state \c a ) (a = 1,4), 
and u Cb represents the probability amplitude of the incoherent continuum state |c&) {b = 5, 8). 
As already explained at the beginning of Sec. II the incoherent continuum states |c&) are 
not presented in Fig. [2] to avoid the complexity of the figure. They have the same quantum 
numbers as the coherent continuum states \c a ), but located at different energies. 

Having solved Eqs. (E])-© and Eqs. (SHE]), we can now calculate the total (angle- 
integrated) ionization yield, R(t), from the relation, 



R(t) = J2Rc a (t), 



(9) 



a=l 



where the partial photoelectron yields R Ca (t), into each coherent and incoherent continuum 
state \c a ), (a = 1,8) can be calculated through the following formulae: 



R ci (t) = f dt'rgKi 2 , 

J — oo 

R C2 (t) = f dt'{r^\ Ul \ 



+ 4 Im 



n? 3 1 1 + — 



r (p > u, i 2 _l r w Ui 



Re(MiMg) + 4 Im 



R a (t)= f A'{ri?>.l 2 + rgH 2 + rgl«4l 

J — oo k 



n% i + 



+ 4 Im 

RcM = f dt'T^\u 2 \ 

J — oo 

R C5 (t) = f dt'T { £ 5 \u 3 \ 

J — oo 

t 



Re(MiO + 4 Im 



nil 1 1 + 



Re(u 2 u* 4 ) 



Rc 6 (t) 
R C7 (t) 



dt' 



— oo 

t 



r w u, I 2 J- r w in 



dt' 



in i z -i- r w in 



3c7 



4c 6 I 



(P) 



4c7 I 



i? C8 (0 = /* dt'rgM 2 . 

J — oo 



(10) 

(11) 

(12) 

(13) 
(14) 
(15) 
(16) 
(17) 



Since the total ionization yield is a sum of ionization into the coherent and incoherent 
continua, it might be rewritten as, 



R(t)= E 



E<(*)+E<(*) 



.a=l 



b=5 



(18) 



=es,ed 



To see the effects of LICS, it is useful to calculate the branching ratio, B, defined as the 
ratio between the partial ionization yield into each coherent continuum ed and es: 

B = j*. (19) 

The total ionization yield given by Eq. fl9]) and the partial ionization yields given by Eqs. 
f fTU|) -f lT7|) are calculated at the end of the pulses. 



B. Photoelectron angular distribution 



For the purpose of calculating photoelectron angular distribution we need equations be- 
fore angle integration. In order to simplify the calculation of the bound-free dipole matrix 
elements we use a partial wave expansion for the continuum of an alkali-metal atom in a 
coupled \{l' s')i'rriji) basis: 

/ 

|k;nv>= Y a llmil (-iy'- 1/2+m " +m ^j2f + l 



V 1/2 f 



\\L;(l's')j'm f } , (20) 



y my m s i —rriji 



l',m t ,,j' 

where k represents the wave vector of photoelectron, ay mi = Ani 1 ' e~ l6 i'Yit m (0, $), and 6i> 
is the phase shift which is a sum of the Coulomb phase shift and the scattering phase shift; 
recall that the prime index indicates a quantum number for the continuum state. 

We are interested in PAD as a function of polarization angle 9 P . If the final spin state of 
the photoelectron is not detected, we have to incoherently sum over the final spin projection, 
m s i. The partial photoelectron yield into a solid angle, f^, defined by the polar angle 0, 
and the azimuthal angle $, can be written as, 



dR(e,$) 



dtdVL^ 



= 0.589tt Y 

mj =±l/2 m s ,=±l/2 



+ 



3=1 
1 



j=3 



3=3 



(0,$) u 3 



(21) 



where 0.5897T is a conversion factor for the appropriate normalization, and the formula is 
valid for both 4pi/ 2 — 6pi/ 2 and 4p 3 / 2 — 6p 3 / 2 systems, for m 3 - = ±1/2 of the initial state. 
After the considerable angular momentum algebra, we obtain expressions for the differential 
ionization widths r[ a ' ma '\Q, $) from states |1) and |3) for the 4p x / 2 — 6pi/ 2 system: 



,(q,+1/2) 



(0,$) 



la , (22) 



(e,$) 



v/2 



V Z E>( a ) P «W /o *\ (a) VZ pMAv c5^( a ) 

-ijSje%(e,$)ei a) -^-R (a) A v ^ *^ (o) 



R^Y 21 (e,<t>)el 



(23) 



and for the differential ionization widths Tj' m (Q, $), from states |2) and |4): 



p(«,+l/2) 

i 



(e,$) 



v/2 

~3~^ 



V 72 " „(a) 



(e,$) 



(24) 



-i?ge^^ (e,$)e;r ) - ^ge^y^e,^;; 



(a) 



i$]e^y 21 (e,$K 



(a) 



^Se^F 2 _ 1 (e,$)e (a 1 ) 



(25) 



where R^l and represent the radial bound-free matrix elements from state \ j) ( j = 1, 4) 
to the continua es and ed, respectively, by laser a (a = p or d), evaluated in atomic units. 
Here with q = 0, ±1, are the spherical components of the polarization vector of laser a, 
namely e^f 1 = cos6 a , and e±* x = =Fsin#Q,/\/2. The laser intensities, I p and Id, are expressed 
in W/cm 2 . For the coherent continuum, the relevant phase shifts are S s = 1.937 and 5d = 



-6.574, which are the sum of the Coulomb phase shifts, 5 



c 



-4.924 and 5% 



-7.551, 



and the scattering phase shifts, 7r/i s = 6.861 and 7ifi d = 0.977 with fii (I = s,d) being 
the quantum defects estimated from the linear extrapolation of the bound Rydberg s and 
d series of the K atom to the continuum energy of interest. Eq. (I2~T]) together with Eqs. 
( I22p - (l25j) gives PAD for the 4pi/ 2 — 6pi/ 2 system with appropriate normalization, so that the 
angle-integrated quantity becomes identical to the total ionization yield calculated with Eq. 

©• 

Similarly the differential ionization widths from states |1) and |3) for the 4p 3 / 2 — 6^3/ 2 
system, are given by, 



r K+i/2) 

3 



(e,$) 



^d(«)Av (p. <fA P ( a ) 2v ^ p( a )Av (p. (j)\p( a ) 
—K jes e r o{U,^)e - -^-j=K jed e r 20 (U,<Pje 



^R^Y 2l (Q^)e[ a) - ^R^Y^Q^l 
15 vl5 



(26) 
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i(a,-l/2) 



(6,$) 



^ R^Y 22 (Q, $)eS a) - -^=R^Y 21 (Q, 

• x 15 



15 



and from states |2) and |4) they are derived as, 



r . 



(",+1/2) 



(©,<&) 



-^£ ) e^loo(e,$)eS Q) - ^=i?ge^r 20 (e,$)eS Q) 



r (a.-l/2) 



(e,$) 



^i?t )el ^oo(e,$)e^ - ^<]e^F 2 o(e,$)e^ 



(27) 



/a, (28) 



(29) 



'15 J V15 

The Eqs. fl22|) - (!29|) are applicable for both probe and dressing lasers. Recalling that the 
polarization vector of the dressing laser is parallel to the quantization axis (9d = 0°), the 
relative polarization angle between the probe and dressing lasers becomes identical to 9 p . 
Because of the symmetry properties of spherical harmonics, we can show that r'^'" v ' ) and 



r 



(a,m s ,) 



is equal to V 



(a-m s ,) 



and T 



(a-m s ,) 



respectively, by interchanging ej°' and e 



(a) 

i ■ 



III. NUMERICAL RESULTS AND DISCUSSION 



In this section we present numerical results and discussions. All the necessary single and 
effective two-photon dipole matrix elements needed for our schemes have been obtained using 
quantum defect theory and Green's function technique. The calculated atomic parameters 
such as Rabi frequencies, asymmetry parameters, ac Stark shifts, for the ipy 2 — 6P1/2 an d 
4j?3/2 — 6p 3 / 2 systems, are listed in Tables I and II, respectively. 

For the K 4pi/ 2 — Gpi/2 system only Rabi frequencies are 6 P dependent, ionization widths 
and ac Stark shifts by the probe laser do not present any 9 P dependence since the initial 
state is isotropic 12|, |l4| (recall that by using a linearly polarized auxiliary laser the ground 
state, 4si/2, of the K atom is excited to the 4pi/2 state, thus all the magnetic sublevels for 
the 4p!/ 2 state are occupied with the same probability). In contrast, for the K 4p 3 / 2 — 6p 3 / 2 
system the atomic parameters such as Rabi frequencies, and ionization widths and ac Stark 
shifts by the probe laser depend on the polarization angle 8 P . This is due to the fact that not 
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all the magnetic sublevels for the initial state, 4p 3 / 2 , are occupied with the same probability. 
The total and partial asymmetry parameters are independent on the laser fields 3i], and, 
obviously, do not depend on the polarization angle of the probe laser for both K 4p x /2 — 6P1/2 
and 4p 3 / 2 — 6p 3 / 2 systems. 

When we solve the set of amplitude equations a special care has to be taken: Note 
that the probability amplitudes of the initially occupied states, |1) and |2), have arbitrary 
phases. However, if we use the amplitude equations, once the initial conditions are given 
for the probability amplitudes u%{t = — oo), % — 1,4, initial coherence inevitably exists due 
to nonzero value of UiU* (with i ^ j). We have to avoid any coherent interference between 
the ionization paths starting from |1) and |2). Therefore, we should separately solve the 
set of amplitude equations Eqs. ©-([S]) with either u\{t = — oo) = 1 and Ui(t = — oo) = 
(for i = 2,3,4), or w 2 (t = — oo) = 1 and Ui(t = — oo) = (for i = 1,3,4), and average the 
photoelectron angular distribution given by Eq. ( 12~T|) over rrij of the initial state: 



dR{Q,$) _ 1 
dtain^ ~ 2 



dtdVt\ 



dR(e,§) 



-1/2 



dtdil-i 



'1/2. 



(30) 



Pulse durations and peak laser intensities are chosen to be r v = 1 ns (FWHM) and I p = 1 
MW/cm 2 for the probe laser, and 10 ns < r d < 15 ns (FWHM) and 100 MW/cm 2 < I d < 500 
MW/ cm 2 for the dressing laser. If the probe and the dres sing pulse durations are comparable 
the LICS resonance profile is going to be smeared out [ll| due to the ac Stark shifts. In 
order to circumvent this problem the pulse duration of the dressing laser was chosen to be 
much longer than that of the probe laser, since, under the condition that r d ^> r p , atomic 
states are quasi-statistically Stark-shifted by the strong dressing pulse during the interaction 
with the probe pulse. By substituting the atomic parameters listed in Tables I and II into 
Eqs. ([2])-(j5]), we can easily solve those equations for the given peak intensities, detunings, 
and temporal profile of the lasers. Once the solution is obtained for Ui(t) (i = 1,4), the 
total and partial ionization yields can be calculated from Eq. (J9]) and Eqs. (|T0|) - (|T7|) . The 
radiative lifetimes of 4p ly / 2 and 4p 3 / 2 , and 6p±/2 and 6p 3 / 2 levels, are about 26 ns and 345 ns, 
respectively, and are included in the numerical calculations. 

In order to check the consistency of our results an alternative formalism based on the 
density matrix equations was used to calculate the dynamics of the system. The amplitude 
equations approach has the advantages of dealing with fewer number of differential equations 
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and obtaining more compact formulae for the ionization yield. On the other hand, the density 
matrix equations approach has the advantage of its capability to control the coherence 
through the off-diagonal density matrix elements (that are proportional to p^- = UiU*, with 
i 7^ j), and therefore it is suitably used for mixed states (such as 4p 1 / 2 (mj = ±1/2) or 
Ap 3 / 2 (mj = ±1/2)) when at least two levels with arbitrary phase are initially occupied. 
Details about the density matrix approach are given in Appendix |Aj and the numerical 
results are, of course, identical to the ones obtained by using the amplitude equations. In 
the following subsections we present numerical results for the K ^p\/2 — Gpi/2 an d 4p 3 / 2 — 6p 3 / 2 
systems. 

A. K Api/2 — Gpi/2 system 

Figs. |3]^a) and 0(b) show the variation of the total ionization yield and branching ratio 
as a function of two-photon detuning, 6, at four different values of the polarization angle, 
6 p = 0°, 30°, 60° and 90°. Note that the polarization angle of the dressing laser is fixed 
to 6d = 0°. Pulse durations and peak laser intensities are chosen to be t p = 1 ns and 
I p — 1 MW/cm 2 , and = 10 ns and Id = 100 MW/cm 2 , for the probe and dressing lasers, 
respectively. Clearly, the profile of the LICS resonance and the branching ratio as a function 
of detuning for the Ap\n — Qpi/2 system changes by varying the polarization angle. The 
position of the LICS resonance is also altered by varying 9 P , specifically the maximum of 
the LICS profile shifts toward larger values of the detuning as 9 P increases and its minimum 
vanishes completely when 9 P = 90°. At 9 P = 90° the value of the asymmetry parameter, 
which is connected to the resonance profile, q = —6.59, is much larger compared to the case 
when both lasers are linearly polarized in the same direction, and the asymmetry parameter 
takes the value, q = —0.91. That particular value of the asymmetry parameter for 9 P = 90° 
is due to the fact that the corresponding angular coefficients for the s and d ionization 
channels are equal, and the radial matrix elements have opposite signs. 

Figs. Ufa) and \MJd) show the variation of the total ionization yield and branching ratio 
as a function of detuning at three different dressing laser intensities Id = 100, 200, and 
500 MW/cm 2 , with the probe laser intensity and the pulse durations fixed to be I p = 100 
MW/cm 2 , t p = 1 ns, and = 15 ns. The polarization angle is 9 p = 30°. As we have already 
seen in Figs. 131(a) and E)(b), the ionization yields and branching ratios vary significantly 
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near resonance. The LICS structure is naturally broadened as the dressing laser intensity is 
increased. 

For the particular polarization geometry shown in Fig. [T]the azimuthal angle dependence, 
$, of the photoelectron signal, through the spherical harmonics Y/ m ($, G), does not vanish 
as it happens when the polarization axes of both lasers are parallel to each other and the 
cylindrical symmetry of PAD is broken. This is due to the presence of the spherical harmonics 
with m ^ in the differential ionization widths formulae Eqs. (T22]) - (T29|) . Before studying 
PAD for LICS it would be instructive to give an answer to the following question: What 
is the modification of PAD for one-photon ionization from the initial state 4px/ 2 through 
the variation of the polarization angle of the probe laser without the dressing laser, i.e., 
Id = 0? Since the initial state, ^p\/2, is spherically symmetric (recall that both rrij = ±1/2 
sublevels are equally populated), one could intuitively guess that the magnitude of PAD does 
not change, and PAD just aligns along the polarization axis of the probe laser. Under the 
condition of Id = 0, three-dimensional (3D) PAD is plotted in Figs. E]^a)-[5|c) as a function 
of photoelectron angles O and $, for the three different values of the polarization angle, 
dp = 0°, 45°, and 90°. As expected, PAD changes its orientation along the polarization 
direction of the probe laser. The fourfold rotational symmetry ($ — > ir + $, and $ — > — $, 
at = 90°) that exists when both lasers are linearly polarized along the quantization axis, 
jreaks into a twofold symmetry ($ — > n + $) when the polarization of the probe laser varies 
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16|. 



Now we consider the case of LICS, i.e., the dressing laser is turned on. 3D PAD's of the 
K 4pi/ 2 — Qpi/2 system at 6 P = 60° are shown in Fig. Et^-IH^c) for the three representative 
detunings corresponding to the far-off resonance (S = —4 GHz), maximum (6 = —0.62 GHz), 
and minimum (5 = 0.44 GHz) of the branching ratio (see Fig. [3]^b)). The view point of all 
3D plots in this paper is from the xy-plane with the Cartesian coordinates (2,2,0), if not 
otherwise stated. At far-off resonance (Fig. |6]^a)), the 3D PAD again tends to follow the 
change of the polarization angle, 8 P = 60°, with some small distortion due to the dressing 
laser. The distortion, however, is almost invisible, since the interference effect through LICS 
is negligible at far-off resonance. In Figs. Mb) and |6]^c), we see that the 3D PAD's are 
significantly modified. Especially in Fig. [6]^c) , a maximum distortion is observed in PAD 
due to the strong destructive interference between the es and ed partial waves. 
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B. K Ap 3/2 - 6p 3 / 2 system 



We now turn to the case of the 4p 3 / 2 — 6^3/2 system. This system is somehow different 
from the 4pi/2-6pi/2 system, because the initial state, 4p 3 / 2 , is not spherically symmetric: 
Only rrij = ±1/2 out of all possible rrij = ±1/2, ±3/2 magnetic sublevels are equally 
occupied by the auxiliary laser, and for this reason a different behavior is expected. 

In Fig. [7(a) we plot the variation of the total ionization yield as a function of two-photon 
detuning at four different values of the polarization angle, 9 P = 0°, 30°, 60°, and 90°. The 
LICS structure in Fig. [7(a) is not quite similar to that plotted in Fig. [3(a), because, 
although gig and q\f are the same for both systems, 513 itself is different. More interestingly, 
the variation of the branching ratios shown in Fig. [71(b) is substantially larger than that 
shown in Fig. [3(b) as 6 P increases. For 9 P = 90° the ionization yield into the ed continuum 
at the detuning close to 8 = —0.9 GHz is almost 40 times enhanced compared to that 
into the es continuum. This suggests that an appropriate choice of the probe polarization 
angle and the two-photon detuning leads to the control of ionization into different channels. 



Recent experiments 



171 ] performed for ionization from an excited state of Xe with linearly 



and circularly polarized lasers have demonstrated that the ionization products into different 
continua can be separated by varying the polarization of lasers. The variation of the total 
ionization yields and the branching ratios as a function of detuning 5 at three different 
dressing laser intensities, Id = 100, 200, and 500 MW/cm 2 , is presented in Figs. [8(a) and 
[8(b) with the rest of the parameters being the same as those in Fig. [H 

In Figs. E(a)-|n(c) we plot the 3D PAD for one-photon ionization from the 4p 3 / 2 state by 
the probe laser at 9 P = 0°, 45°, and 90°, without the dressing laser, i.e., Id = 0. Compared 
to the 4pi/2 — Gpi/2 system (see Figs. E(a)-[5(c)), PAD's drastically change the shape with 
the detuning when the probe polarization angle is varied. 

Now we return to the case for LICS by turning on the dressing laser, and present the 3D 
PAD's for the 4p3/ 2 — 6p 3 / 2 system in Fig. fIDT a)- fT0T c). at 9 P = 60°, for three representative 
detunings corresponding to the far-off resonance (5 = —4 GHz), maximum (5 = —0.89 GHz), 
and minimum (5 = 0.42 GHz) of the branching ratio (see Fig. [7(b)). The modification of 
the 3D PAD's, presented in Figs. [9] and [10], is more than we expect: The variation of the 
sidelobes of 3D PAD's at different polarization angles, which are absent for the 4pi/ 2 — Qpi/2 
system, is striking. The sidelobes are due to the ionization into the continuum (that 
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is inaccessible through one-photon ionization from the 4pi/2 state), and are present even if 

e p = o°. 

As we have already noticed, not only the photoelectron angular distribution but also the 
angle-integrated ionization yield is affected by the relative polarization angle between the 
probe and dressing lasers. This effect is called linear dichroism and is very attractive from 
the experimental point of view, since it is much easier to measure the total ionization yield 
than the PAD. Linear dichroism can be experimentally used to determine the ratio of the 
dipole matrix elements into the different continua or the relative phase shift between 
the partial waves of the continua. The normalized linear dichroism is defined as I9I ]. 

rn _ R(0 P = 90°) - R(6 P = 0°) 

R(9 P = 90°) + R(9 P = 0°)' { ' 

where R{9 p = 0°) and R{6 P = 90°) represent the total ionization yield, or equivalently 
the angle-integrated photoelectron signal when the polarization axis of the probe laser is 
parallel and perpendicular with respect to that of the dressing laser, respectively. In Figs. 
flTT a) and fTTT b) we plot the linear dichroism, LD, as a function of two-photon detuning for 
the K 4pi/2 — 6P1/2 and 4p 3 / 2 — 6p 3 / 2 systems, respectively. The magnitude of the linear 
dichroism changes drastically for both systems around the LICS resonance, and it shows a 
large maximum for the two-photon detunings around the deep LICS minimum at 9 P = 0°. 



IV. SUMMARY 



In this paper we have theoretically investigated the effects of the relative polarization 
angle between the probe and dressing lasers on the total (angle-integrated) ionization yield, 
branching ratio, and PAD through LICS for the K 4pi/ 2 — Qpi/2 arid 4p 3 / 2 — 6p 3 / 2 systems 
in a particular geometry with both probe and dressing lasers being linearly polarized. Am- 
plitude equations and alternatively density matrix equations formalisms have been used to 
study the dynamics of the ionization process. We have shown that the ionization yield and 
the branching ratio are strongly dependent on the relative polarization angle between the 
lasers. Moreover we have found that ionization into the different continua, branching ratios, 
and PAD are significantly altered by the change of the polarization angle. Our findings sug- 
gest that the relative polarization angle can be another doorknob to control the ionization 
dynamics through LICS. We have also calculated linear dichroism for the angle- integrated 
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ionization yield, which turned out to be quite large at the two-photon detunings close to the 
LICS minimum. 



Acknowledgments 

GB acknowledges financial support from Japan Society for the Promotion of Science 
(JSPS). The work by TN was supported by the Grant-in- Aid for scientific research from the 
Ministry of Education and Science of Japan. 



Appendix A: Time-dependent density matrix equations 



Based on the density matrix approach 20(, we study the temporal evolution of the K 
atom in the laser field given by Eq. (JTJ) . Briefly we solve the following set of time-dependent 
differential equations for the slowly varying density matrix a: 



&u = -Ticra - 21m 
—Tj<jjj + 21m 



a 



a 



j.i 



j=3 \ 
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i=l 



IjiJ 

i 
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„• \ 4 / • \ 



'13 
1 



i'=l 



(Al) 
(A2) 
a fj ,(A3) 



= — + rv)^ + i h — o-^ — ^^iv, i ay , 



3=3 
2 



9ii 



3=3 
2 



E% 1- 



8=1 



(A4) 
(A5) 



where the indices take the following values i, %' = 1, 2 and = 3, 4, with j ^ j' and j ^ f. 
All the density matrix elements for the continuum have been adiabatically eliminated from 
the Eqs. ( 1All) -( fA5l) . Note that we have used the rotating wave approximation and the 
slowly varying density matrix elements to derive the above equations: an = pa, (i = 1,4), 
(Tij = Pije- iSstati ->v\ (i = 1,2, and j = 3,4), a ic = p ic e~ iuJat , [a = p or d, and % = 1,4), 
where Pij(t) = Ui(t)u*(t) are the density matrix elements. 5ij is the two-photon detuning 



defined by 5. 



1.1 



5 



static, ij 



+ 5, 



stark, ij j 



where the static detuning is defined by 5. 



static, ij 



(Ei + Tiujp) — (Ej + hud), and 5 sta rk, ij is the total dynamic ac Stark shift defined by S, 



stark, ij 



(Si + Sl a> ) — (Sf + Sj a> ). Now the above set of density matrix equations is solved with 



7 (p) 



7 (<9 
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the following initial conditions: au(t = — oo) = 1/2, and <Jjj(t = — oo) = <Jij(t = — oo) = 0, 
for i = 1,2 and j = 3,4. The total (angle-integrated) ionization yield is derived as, 

2 4 2 4 



= / dt> £ r?V« + £(rf + if K, + 4 £ £ im 

J - 00 li=l j=3 i=lj=3 



In, 



Re(<Jij) 



(A6) 



The probability that a photoelectron is ejected into a solid angle, fl-^, is given by the following 
formula: 



dR(Q,$) 



a 



= 0.589 vr £ £rS w) (e,$)a, 

m s /=±l/2 U=l 

+ £ [rf "^(e, $) + rj w) (e, $) 
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2 4 

+ 2Re 

+ 2Re 
+ 2Re 
+ 2Re 



E E Vr! w He, $) V(ri w) (e, $)) V, 

i=l j=3 





) (e,$))*< 







r 3 ^' ) (e,$)J(r 4 ^' ) (e,$)) a 34 



(A7) 



which can be shown to be equivalent with Eq. ( I30p . The numerical results for both K 
4pi/2 _ Gpi/2 and 4p 3 / 2 — 6^3/2 systems obtained in the density matrix and amplitude equations 
formalisms are, of course, identical. 



Appendix B: Ionization widths 

The partial ionization width from state \j) to the continuum |c) produced by laser a, is 
defined as: 

rg' ms,) (0,$) = 2Tt\Df/ ns 'X®^)\\ (Bl) 

where D^'™"'^ = —E a (t) S q =±i,o( c l r g e g a ' > b) = /jj"' ms '' represents the one-photon dipole 
matrix element between states \j) and |c), expressed in the length gauge, and calculated at 
energy E c = Ej + u a . The total ionization width integrated over the solid angle, f^, defined 
by the polar angles (9, <3>), of the ejected photoelectron is given by, 

if } = E Er£' ms,) - (B2) 

m„/=±l/2 c 
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The summation over c implies that the summation is taken over the all allowed continuum 
states. 



Appendix C: Ac Stark shifts 



The dynamic ac Stark shift of the energy of state \j) due to both bound and free states 
\k), caused by laser a, is, 



-,(<*) 



E E 

m s /=±l/2 k 



\D 



0>"V)|2 
jk 



+ 



\D 



0>«V)|2 
jk 



Ej + oj a — Ek + is Ej — u> a — Ek + is 



(CI) 



The sum here contains both summation over the bound and integration over the continuum 
states, and e is an infinitely small number. 



Appendix D: Two-photon Rabi frequency 

The total two-photon Rabi frequency Qij between the state \i) and \j) is given by, 

where the partial two-photon Rabi frequency between state \i) and \j) coupled through the 
continuum |c) with the energy E\ + u p ~ £3 + uia is defined as, 

Q. c , ( 1 - — ) = V / d£ c -^ . (D2) 

The imaginary part of the partial Rabi frequency is connected to the partial asymmetry 
parameter qfj, and is given by, 



Qij m s /=±l/2 



(D3) 
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Table I. Atomic parameters for the K 4P 1 / 2 -6pi/ 2 system. Q is measured in 
rad/s, T in s -1 , S in rad/s, and Id in W/cm 2 . 
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Table II. Atomic parameters for the K 4P3/ 2 -6£>3/ 2 system. Q is measured in 
rad/s, T in s -1 , S in rad/s, and /, in W/cm 2 . 
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FIG. 1: Quantization axis and the polarization vectors, and e^ d \ for the probe and dressing 
lasers defined for this work. The polarization vector lies in the xz-plane, and the quantization 
axis is taken along the z axis. 
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FIG. 2: Color online. Level scheme considered in this paper for the K 4p]y 2 — &Pi/2 system. 
Depending on the polarization angle of the probe laser, 8 p , different transition paths have to be 
considered, (a) 9 P = 0°, (b) 9 P = 90°, and (c) 0° < 8 P < 90°. Similar level scheme can be drawn 
for the K 4p 3/ / 2 — 6p 3 / 2 system. 
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FIG. 3: Color online, (a) Total ionization yield and (b) the branching ratio between the partial 
ionization yields into each es and ed continuum for the K 4p 1( / 2 — €>Pi/2 system as a function of 
two-photon detuning 5. Pulse durations and peak laser intensities are chosen to be r p = 1 ns 
and I p = 1 MW/cm 2 , and = 10 ns and 1^ = 100 MW/cm 2 , for the probe and dressing lasers, 
respectively. The polarization angle takes the values of 6 P = 0°,30°,60° and 90°. 
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FIG. 4: Color online, (a) Total ionization yield and (b) the branching ratio between the partial 
ionization yields into each es and ed continuum for the K 4p 1( / 2 — ®Pi/2 system as a function of 
two-photon detuning 5 for the three different dressing laser intensities, Id = 100, 200, and 500 
MW/cm 2 . The intensity of the probe laser is Id = 1 MW/cm 2 . Pulse durations are chosen to be 
t p = 1 ns and = 15 ns, for the probe and dressing lasers, respectively. The polarization angle is 
6L = 30°. 
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FIG. 5: Three-dimensional photoelectron angular distribution due to one-photon ionization from 
the K Ap l /2 state by the probe laser field only, at three different polarization angles 9 P = 0°, 45°, 
and 90°. Pulse duration and peak intensity are t p = 1 ns and I p = 1 MW/cm 2 for the probe laser. 
The view point is from the positive y-axis. 
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(a) 




FIG. 6: Three-dimensional photoelectron angular distribution for the K 4p ly / 2 — Qpi/2 system at 
three different two-photon detunings 5 = —4, —0.62, and 0.44 GHz. Pulse durations and peak 
intensities are t p = 1 ns and I p = 1 MW/cm 2 for the probe laser, and = 10 ns and I p = 100 
MW/cm 2 for the dressing laser. The polarization angle is 6 P = 60°. 
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FIG. 7: Color online, (a) Total ionization yield and (b) the branching ratio between the partial 
ionization yields into each es and ed continuum for the K 4p 3 / 2 — 6^3/2 system as a function of 
two-photon detuning 5. Pulse durations and peak laser intensities are chosen to be r p = 1 ns 
and I p = lMW/cm 2 , and = 10 ns and Id = 100 MW/cm 2 , for the probe and dressing lasers, 
respectively. The polarization angle takes the values 9 P = 0°,30°,60° and 90°. 
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FIG. 9: Same as in Fig. [5] but for the K 4j> 3 /2 — 6j> 3 / 2 system. The view point is from the positive 
y-axis. 
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FIG. 10: Three-dimensional photoelectron angular distribution for the K 4p 3 / 2 — 6^3/2 system at 
three different two-photon detunings 5 = —4, —0.89, and 0.42 GHz. Pulse durations and peak 
intensities are t p = 1 ns and I p = 1 MW/cm 2 for the probe laser, and = 10 ns and I p = 100 
MW/cm 2 for the dressing laser. The polarization angle is 6 P = 60°. 
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FIG. 11: Linear dichroism, LD, for the (a) K 4p 1 / 2 — Qpi/2 system and (b) K 4p 3 / 2 — 6p 3 / 2 system 
as a function of two-photon detuning 5. Pulse durations and peak intensities are t p = 1 ns and 
I p = 1 MW/cm 2 for the probe laser, and = 10 ns and I p = 100 MW/cm 2 for the dressing laser. 
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